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ABSTRACT 

We propose a scenario in which massive stars form in a self-gravitating gaseous disc 
around a supermassive black hole. We analyze the dynamics of a disc forming around a 
supermassive black hole, in which the angular momentum is transported by turbulence 
induced by the disc's self-gravity. We find that once the surface density of the disc 
exceeds a critical value, the disc fragments into dense clumps. We argue that the 
clumps accrete material from the remaining disc and merge into larger clumps; the 
upper mass of a merged clump is a few tens to a few hundreds of solar mass. 

This picture fits well with the observed young stellar discs near the SgrA* black 
hole in the Galactic Center. In particular, we show how the masses and spatial dis- 
tribution of the young stars, and the total mass in the Galactic Center discs can be 
explained. However, explaining the origin of the several young stars closest to the black 
hole (the S-stars) is more problematic: their orbits are compact, eccentric, and have 
random orientation. We propose that the S-stars were born in a previous starburst(s), 
and then migrated through their parent disc via type I or runaway migration. Their 
orbits were then randomized by the Rauch-Tremaine resonant relaxation. 

We then explore the consequences of the star-formation scenario for AGN discs, 
which are continuously resupplied with gas. We argue that some compact remnants 
generated by the starburst will get embedded in the disc. The disc-born stellar-mass 
black holes will interact gravitationally with the massive accretion disc and be dragged 
towards the central black hole. Merger of a disc-born black hole with the central black 
hole will produce a burst of gravitational waves. If the central black hole is accreting at 
a rate comparable to the Eddington limit, the gas drag from the accretion disc will not 
alter significantly the dynamics of the final year of merger, and the gravitational waves 
should be observable by LISA. For a reasonable range of parameters such mergers 
will be detected monthly, and that the gravitational-wave signal from these mergers 
is distinct from that of other merger scenarios. Also, for some plausible black hole 
masses and accretion rates, the burst of gravitational waves should be accompanied 
by a detectable change in the optical luminosity of the central engine. 



1 INTRODUCTION 



It is widely believed that self-gravitating accretion discs can 
form around supermassive black holes (SBH) in AGNs. The- 
oretical models show that the AGN accretion discs must 
become self-gravitating if they extend beyond a fraction 
of a parsec away from a central black hole (Paczynski 
1978, Kolykhalov and Sunyaev 1980, Schlosman and Begel- 
man 1987, Kumar 1999, Jure 1999, Goodman 2003). Self- 
gravitating discs are unstable to fragmentation on a dynam- 
ical timescale; self-gravity of AGN accretion discs is a major 
issue in understanding how gas is delivered to the central 
black hole. It is likely that star formation will occur in the 
outer parts of an AGN accretion disc (Kolykhalov and Sun- 
yaev 1980, Schlosman and Begelman 1987). 

There are two lines of observational evidence that a 



starburst within SBH's radius of influence may be a common 
phenomenon: 

(i) Levin and Beloborodov (LB, 2003), using pre- 
existing data, have identified a disc of young massive stars 
which are moving clockwise in the gravitational potential 
of SgrA*, the SBH at the center of the Galaxy. They have 
argued that the stars were born in a dense accretion disc 
which existed several million years ago in the Galactic Cen- 
ter. The presence of the disc was later confirmed by Genzel 
et. al. (2003) and Paumard et al. (2005), who have used up- 
dated data sets; these authors also found a second counter- 
clockwise disc of young massive stars. For the past few years, 
an alternative possibility for the origin of the young stars 
was considered on equal footing with LB'S proposal. In the 
alternative scenario (Gerhard 2001, Hansen and Milosavlje- 
vic 2003), the stars were born tens of parsecs from SgrA* 
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and were originally members of a massive star cluster. This 
cluster then spiraled in towards SgrA* due to dynamical 
friction with an inner bulge, and deposited the stars near 
SgrA*. However, this scenario is now disfavored due to two 
recent observations: (a) In the cluster scenario one expects 
many young low-mass stars to also be present in the Galac- 
tic Center and for them to produce copious amounts of x- 
ray emission from their coronae. However, Chandra obser- 
vations show that the Galactic Center is a relatively weak 
x-ray emitter, and thus the presence of a multitude of low- 
mass stars is ruled out (Nayakshin and Sunyaev, 2005); (b) 
Paumard et al. (2005) find no early-type stars outside of 
the central half-parsec, which again strongly argues against 
the in-spiraling cluster scenario. Therefore, most likely the 
young stellar disc formed in situ as result of fragmentation 
of a gaseous disc, as was argued by LB. 

(ii) Bender et al. (2005) have used HST to identify a 
compact disc of A-stars which is located deep in the gravi- 
tational potential of a SBH in the nucleus of M31. The stars 
are likely much younger than the SBH (this is true unless 
one assumes that the SBH is older than ~ 10 s years), and 
the strong tidal barrier makes it unlikely that the disc is a 
remnant of a tidally disrupted star cluster since any clus- 
ter would get disrupted at a larger distance from the SBH 
(Nayakshin, 2005). Therefore, the disc of A-stars is most 
likely the remnant of a gaseous accretion disc which existed 
in M31 about a hundred million years ago. 

Star formation in an SBH radius of influence may 
be connected to supporting the high accretion of some 
AGNs (Thompson et al. 2005) and at the same time may 
help explain low luminocity of others (Tan and Blackman 
2005). The dynamics of the fragmenting disc is strongly 
affected by the feedback energy input from the starburst. 
Collin and Zahn (1999) have conjectured that the feed- 
back from this star formation will prevent the accretion disc 
from becoming strongly self-gravitating. However, Goodman 
(2003) has used general energy arguments to show that the 
feedback from star formation is insufficient to prevent an 
AGN disc with the near-Eddington accretion rate from be- 
coming strongly self-gravitating at a distance of 10 4 — 10 5 
Scwartzschild radii from the central black hole (about O.lpc 
for our Galactic Center). This is distinct from the case 
of galactic gas discs, for which there is evidence that the 
feedback from star formation protects them from their self- 
gravity. 

In this paper we concentrate on the physics of the self- 
gravitating disc and make a semi-analytical estimate of the 
possible mass range of stars formed in such discs (Sections 
2 and 3). Our principal conclusion is that the stars can be 
very massive, up to hundreds of solar masses. This conclu- 
sion is in qualitative agreement with two recent independent 
observations: Nayakshin and Sunyaev (2005) and Paumard 
et al. (2005) have shown that the young stars in the GC 
must were produced in a starburst with the top-heavy IMF 
strongly favoring massive stars. In Section 4 we specialize to 
the case of SgrA* discs and show that the mass and the col- 
umn density distribution of the marginally fragmenting disc 
are consistent with those of the currently observed stellar 
discs. We also address the puzzle of the several the young 
stars in the central arc-second (the S-stars). Their orbits 
present a problem for the disc-starburst picture. Their ex- 
treme proximity to SgrA*, eccentric orbits, and random in- 



clinations exclude the possibility that they were born from 
the gaseous disc at their current location. We argue that 
instead they were born in a disc at a larger distance from 
SgrA*, but then migrated inwards due to gravitational in- 
teraction with the disc. Their eccentricity and inclination 
angles were randomized by relatively fast resonant relax- 
ation, a process discovered by Rauch and Tremaine in 1996. 

In section 5, we consider a self-gravitating AGN accre- 
tion disc which is continuously supplied by gas on a timescale 
greater than the lifetime of massive stars. Some of the black- 
hole remnants of the stars become embedded in the disc and 
and migrate inward on the timescale of ~ 10 7 years. The 
merger of the migrating black hole with the central black 
hole will produce gravitational waves. We show that for a 
broad range of AGN accretion rates the final inspiral is unaf- 
fected by gas drag, and therefore the gravitational-wave sig- 
nal should be detectable by LISA. The rate of these mergers 
is uncertain, but if even a fraction of a percent of the disc 
mass is converted into black holes which later merge with 
the central black hole, then LISA should detect monthly a 
signal from such a merger. The final inspiral may occur close 
to the equatorial plane of the central supermassive hole and 
is likely to follow a quasi-circular orbit, which would make 
the gravitational-wave signal distinct from those in other as- 
trophysical merger scenarios. If the disc-born black hole is 
sufficiently massive, it will disrupt accretion flow in the disc 
during the final year of its inspiral, thus making an optical 
counterpart to the gravitational-wave signal. 



2 PHYSICS OF A FRAGMENTING 
SELF-GRAVITATING DISC 

The importance of self-gravitating accretion discs in astro- 
physics has long been understood (Paczynski 1978, Lin and 
Pringle 1987). It was conjectured that the turbulence gen- 
erated by the gravitational (Toomre) instability may act 
as a source of viscosity in the disc. This viscosity would 
both drive accretion and keep the disc hot; the latter would 
act as a negative feedback for the Toomre instability and 
would keep the disc only marginally unstable. Recently, 
there has been big progress in our understanding of the self- 
gravitating discs, due to a range of new and sophisticated 
numerical simulations (Gammie 2001, Mayer et. al. 2002, 
Rice et. al. 2003). In our analysis, we shall rely extensively 
on these numerical results. 

Consider an accretion disc which is supplied by a gas 
infall. This situation may arise when a merger or some other 
major event in a galaxy delivers gas to the proximity of a 
supermassive black hole residing in the galactic bulge. Let 
E(r) be the surface density of the disc. We follow the evolu- 
tion of the disc as S(r) gradually increases due to the infall. 

We begin by assuming that initially there is no viscosity 
mechanism, like Magneto- Rotational Instability (MRI), to 
transport the angular momentum and keep the disc hot*. 
This assumption is valid when the ionization fraction of the 
gas in the disc is low, i.e. when the gas is far enough from a 

* When the disc begins to fragment, the viscosity due to self- 
gravity-driven shocks exceeds the one due to MRI; see below. 
Therefore, while computing the disc parameters at fragmentation, 
it is reasonable to ignore MRI 
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Figure 1. The temperature of the critically fragmenting disc as 
a function of the orbital period. 



central source (about a thousand Schwarzschild radii from 
the supermassive black hole). We also, for the time being, 
neglect irradiation from the central source; this may be a 
good assumption if a disrupted molecular cloud forms a disc 
but the accretion onto the hole has not yet begun. As will 
be discussed below, irradiation is important for some regions 
in the AGN discs we are considering. However, as shown in 
the following subsection, inclusion of irradiation or other 
source of heating will only strengthen the case for formation 
of massive stars. 

We assume, therefore, that the forming disc cools until 
it becomes self-gravitating; this happens when 



Q 



(1) 



Here c s is the isothermal speed of sound at the mid-plane 
of the disc, and Q is the angular velocity of the disc. Nu- 
merical simulations show that once the disc becomes self- 
gravitating, turbulence and shocks develop; they transport 
angular momentum and provide heating which compensates 
the cooling of the disc (Gammie 2001). We thus assume, in 
agreement with the simulations, that when the disc exists, 
it is marginally self-gravitating, i.e. Q = 1. Then 



ttGS 



and 



T ~ 2m p c a /kB 
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Here T is the temperature in the mid-plane of the disc, m p 
is the proton mass, and fcs is the Boltzmann constant. 

The one-sided flux from the disc surface is given by the 
modified Stephan-Boltzmann law: 



(4) 



where T e g is the effective temperature. It is related to the 
mid-plane temperature by 



T 4 f(r) 



(5) 



where r = kS/2 is the optical depth of the disc; here k(T) 
is the opacity of the disc. The function /(r) = r for an 
optically thin disc, and /(r) = 1/t for an optically thick 
disc. We combine these two cases in our model by taking 

T 



f(r) 



+ 1 



(6) 



We have used Eq. (3) in the last step of Eq. (5). The flux 
from the disc is powered by the accretion energy: 

a(nCr4- (7) 



3 9 • 

F= —Q, 2 M = 
8tt 



Here a parametrizes viscous dissipation due to self-gravity 
(Gammie 2001); we have used Eq. (2) in the last step. Using 
Eqs (4), (7), and (5), we can express the viscosity parameter 
a as 

S 5 



8a / 2m p 
a = — I 



9 V fe f 



(*<?)7(t) 



(8) 



It is very important to emphasize that in this model a is only 
a function of £ and fi: the temperature in the mid-plane is 
determined by Eq. (3), and this temperature sets the opacity 
which in turn determines the optical depth r = kE/2. The 
opacity in the range of temperatures and densities of interest 
to us is set by light scattering off ice grains and, in some 
cases, by scattering off metal dust. The relevant regimes 
are worked out in the literature on protoplanetary discs; we 
use the analytical fit to the opacity (in cm 2 /gm) from the 
appendix of Bell and Lin, 1994. 



k = 0.0002 x Tk 
t = 2 x 10~ 16 TZ 7 



for T < 166K, 

for 166K < T < 202K 

for T > 202K. 



(9) 



In the first interval the opacity is due to ice grains; in the sec- 
ond interval the ice grains evaporate, and the opacity drops 
sharply with the temperature; in the third interval (the high- 
est T) dust grains are the major source of the opacity. 

From Eq. (8) we see that as E of the disc increases due 
to the merger-driven infall, the effective viscosity will reach 
Qcrit ~ 1. At this stage the cooling time of the disc becomes 
comparable to the orbital period. Gammie's simulations 
show that in this case the turbulence induced by self-gravity 
is no longer able to keep the disc together, and the disc frag- 
ments. Gammie's simulations give a cr i t — 0.3; the numerical 
value we quote disagrees with Gammie's, but agrees with 
Qcrit quoted by Goodman (2003) since like Goodman wc 
use isothermal speed of sound for a-prescription. Gammie's 
results, although obtained for razor-thin discs, justify the 
key assumption of our model: the disc exists as a whole for 
ct < a C rit; and fragments once a = a cr it- Similar criterion 
for the disc fragmentation was already used in Shlosman and 
Bcgclman, 1987. 

We shall refer to the mid-plane temperature and surface 
density of the marginally fragmenting disc with a — a CI - lt as 
the critical temperature T cr i t and the critical surface density, 

S cr it . 
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We now find the critical surface density and mid-plane 
temperature as a function of Q. We use Eq. (3) to express E 
as a function of T and fi, then substitute this function into 
Eqs. (6) and (8), and set a — a cr it- After simple algebra, we 
obtain 



ff + = q, 
where 

p = 2 



(10) 



TVG 



q = 



K(Tcrit) / kBT CIl 

32a (ttG) 2 



(mp/fe s ) 2 T c 2 rit . 



9a cr it K(T cr it) 
There is an analytical solution to Eq. (10): 

£7 = w — p/(3w), 
where 

™ = {g/2 + [(g/2) 2 + (p/3) 3 ] 1/2 } 1/3 - 



(11) 

(12) 
(13) 



In Figure 1 we make a plot of T cr i t as a function of the 
orbital period, for concreteness we set a cr it = 0.3. The criti- 
cally self-gravitating disc is optically thin if the second term 
of the LHS of Eq. (10) is dominant, and optically thick oth- 
erwise. This can be expressed as a condition on the critical 
temperature: the disc is optically thin if 



T crit < 12K(a crit /0.3) 2/15 , 



(14) 



and optically thick for higher critical temperatures. The an- 
gular frequency above which the critically unfragmented disc 
becomes optically thick is 



fit 



ransition 



16.3 x 10" sec" 



(15) 



We use Eqs. (2) and (3) to find the critical surface den- 
sity E cr u, which is plotted in Fig. 2, and the scaleheight 
/icrit = c s /Q. of a marginally fragmenting disc. The Toomre 
mass M c \ — E cr itfe cr ; t is the mass scale of the first clumps 
which form in the first stage of fragmentation. In Fig. 3, we 
plot the Toomre mass of the critically fragmenting disc as a 
function of the orbital period. 

The value of M c \ is not large enough for the initial clump 
to open a gap in the accretion disc. The newly-born clump 
will therefore accrete from the disc. The Bondi-Hoyle esti- 
mate of the accretion rate gives M c \ ~ f2M c i, i e. we ex- 
pect the mass of the new clump to grow on the dynamical 
timescale until it becomes large enough to open a gap in 
the gas disc. The upper limit M c \ of this value is the mass 
which opens a gap in the original gas disc with E = E cr i t 
just before it fragments: 



Me ~ Mei [127r(Q crit /0.3)] 1/2 (r/h ciit ) 



a/2. 



(16) 



see Eq. (4) of Lin and Papaloizou (1986). Once the gas is 
depleted from the disc, we expect the initial distribution of 
the clump masses to be concentrated between M c i and M c \. 
The clump masses will evolve when clumps begin to merge 
with each other; this is addressed in section III. 

2.0.0.1 Effect of irradiation and other sources of 
heating . So far in determining the structure of the self- 
gravitating disc, we have neglected external or internal heat- 
ing of the disc. This is certainly a poor approximation in 
many cases. Irradiation from AGN or surrounding stars, or 
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Figure 2. The surface density of the critically fragmenting disc 
as a function of the orbital period. 



feedback from star formation inside the disc can be the dom- 
inant source of heating of the outer parts of AGN discs (eg, 
Shlosman and Begelman 1987, 1989). For example, irradi- 
ation from circumnuclear stars will keep the disc tempera- 
ture at a few tens of Kelvin, which is higher than the critical 
temperature we obtained for a self-gravitating disc beyond 
O.lpc away from 1O 7 M black hole; see our discussion for 
the Galactic Center disc in section 4. 

The extra heating will always work to increase the crit- 
ical temperature at which the disc fragmentation occurs. 
Therefore, the values of the critical surface density E cr , 
scaleheight h„, and the mass of the initial fragment M c i 
obtained above should be treated as lower bounds of what 
might be expected around real AGNs or protostars. Higher 
values of these quantities would only strengthen main con- 
clusions of this paper. LB have found that when the rate 
of accretion M is constant throughout the disc, the Toomre 
mass Mci is given by 



= 1.8M, 



Mc 2 



M 



iodd V 3 x 1O 6 M 
and that the gap-opening mass is 

-0.5 



0.2pc 



(&)" 



x 



M 



Mc 2 



L cdd V 3 x 10 6 M G 



0.2pc 



(17) 



(18) 

/ r \0.5 
\Wl) ' 



Here, M is the black-hole mass, and r and h are the radius 
and the scaleheight of the disc. Thus, even if the clumps 
do not merge with each other and stop their growth at the 
gap-opening mass, for high accretion rate the mass of the 
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formed stars will be biased towards the high-mass end. In 
the next section we discuss the effects of clump mergers and 
the mass range of stars born after the disc fragments. 



3 EVOLUTION OF THE FRAGMENTED DISC 

Gammie's simulations show that once the disc fragments, 
the clumps merge and form significantly larger objects. In 
fact, his razor-thin shearing box turned into a single gas 
lump at the end of his simulation. 

For merger to be possible, the clumps should not col- 
lapse into individual stars before they can coalesce with each 
other. Let's check if this is the case. 

Consider a spherical non-rotating clump of radius R and 
mass Mci. First, assume that the clump is optically thin. The 
energy radiated from the clump per unit time is 

Wcooi ~ aT 4 i? 2 «;(T)E ~ M cl oT A n(T). (19) 

This radiated power cannot exceed the clumps gravitational 
binding energy released in free-fall time, G 16 M 25 R~ 25 . To- 
gether with k(T) oc T 2 [since ice grains dominate opacity for 
the optically-thin marginally fragmenting disc-see Eq. (14)], 
this condition implies that 

T < f = T R~ 5/12 , (20) 

where T is a constant for the collapsing optically thin 
clump. The temperature T in Eq. (28) is less than the virial 
temperature, which scales as RT 1 . Therefore, after the col- 
lapse commences, the clump is not virialized while it is opti- 
cally thin. The temperature cannot be much smaller than T, 
since otherwise the cooling rate would be much smaller than 
the rate of release of the gravitational binding energy, and 
the gas would heat up by quasi-adiabatic compression. The 
inequality in Eq. (28) should be substituted by an approx- 
imate equality, and therefore we have during optically-thin 
collapse 

TcxR~ 5/12 . (21) 
The optical depth scales as 

r oc R- 11/3 . (22) 

and hence rises sharply as the clump's radius decreases; as 
the clump shrinks it becomes optically thick^. It is possible 
to show that once the clump is optically thick, it virializes 
quickly with it's temperature T oc R~ . For k oc T 2 (ice 
grains), the cooling time of an optically thick clump scales 
with the clump radius as 

tcooi oc R~ 3 . (23) 

The characteristic timescale for the clump to collide with 
another clump scales with the clump radius as 

^collision oc R 2 . (24) 

t The contraction of an optically thin clump may be complicated 
by sub-fragmentation, since the Jean's mass for such clump scales 
as R 7 / 8 oc t~ 0,23 . We suspect that in most cases the clump be- 
comes optically thick before it sub-fragments, since the Jean's 
mass has a slow dependence on the optical depth. However, only 
detailed numerical simulations can resolve these issues. 



From Eqs. (23) and (24), we see that the collision rate de- 
creases less steeply than the cooling rate as a function of the 
radius of an optically thick clump. Therefore, merger can be 
an efficient way of increasing the clump's mass. 

This conclusion is no longer valid when the tempera- 
ture of the clump becomes larger than ~ 200K; then the 
opacity is dominated by metal dust with k oc T 1/2 . In this 
case the cooling time scales as i CO oi oc R~ 15 . The collision 
timescale increases faster than the cooling time as the clump 
shrinks, and naively one would expect that mergers may not 
be efficient in growing the clump masses. 

However, we have neglected the rotational support 
within a clump. Each clump is initially rotating with an- 
gular frequency comparable to the clump's inverse dynami- 
cal timescale; for example, in a Keplerian disc each clump's 
initial angular velocity is ~ Q/2. Therefore each clump will 
shrink and collapse into a rotationally supported disc, and 
the size of this disc is comparable to the size of the original 
clump (this picture seems to be in agreement with Gammie's 
simulations). Thus rotational support generally slows down 
the collapse of an individual fragment and makes mergers 
between different fragments to be efficient. 

Magnetic braking is one of the ways for the clump to 
lose its rotational support^ (see, e.g., Spitzer 1978). One 
generally expects a horizontal magnetic field to be present 
in a differentially rotating disc due to the MRI (Balbus and 
Hawley, 1991). Ionization fraction in the disc is expected 
to be small, so the magnetic field is saturated at a sub- 
equipartition value B = /3B cq , with /3 << 1. Horizontal 
magnetic field will couple inner and outer parts of the dif- 
ferentially rotating clump on the Alfven crossing timescale 
iaifvon ~ idynamicai//3, and the collapse will proceed on this 
timescale as well. 

What is the maximum mass that the clump can achieve? 
This issue has been analyzed for the similar situation of a 
protoplanetary core accreting from a disc of planetesimals 
(Rafikov 2001 and references therein). The growing clump 
cannot accrete more mass than is present in it's "feeding 
annulus". This gives the maximum "isolation" mass of a 
clump: 

Mis ~ (27rr ^ )3/2 = 27vV2M cl (r/h clit f 2 , (25) 

where, as above, M c \ — E cr itft 2 rit is the mass scale of the 
first clumps to form from a disc; see, e.g., Eq. (2) of Rafikov 
(2001). However, numerical work of Ida and Makino (1993) 
and analytical calculations of Rafikov (2001) indicate that 
the isolation mass may be hard to reach. The consider a 
massive body moving on a circular orbit through a disc of 
gravitationally interacting particles, and they find that when 
the mass of the body exceeds some critical value, an annular 
gap is opened in the particle disc around the body's orbit. 
We can idealize a disc consisting of fragments as a disc of 
particles of a typical fragment mass Mf T . Once a growing 
clump opens a gap in a disc of gravitationally interacting 
fragments, the clump's growth may become quenched. This 
gap-opening mass of the clump M gap is given by Eq. (25) of 
Rafikov (2001): 



t Another way is via collisions with other clumps. 
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Figure 3. The Toomre mass of the critically fragmenting disc as 
a function of the orbital period. 



M is 2 7 /6 7r l/2 \z r 2) { Er 2) 



1/2 



(26) 



We use the numerical factor 2 7 ^ 6 7r l / 2 I = 1.5 appro- 
priate for thin discs. By taking Q = 1 we get 



M gap ~ UM ir (M b /M c i) 1/3 (r/h clit ). 



In Figure 4 the masses M- ls and M gap are plotted as a func- 
tion of radius for a 3 x 1O 6 M0 black hole; when we calculate 
-Mgap we conservatively set Mf r = M c \ and not to the larger 
value M c i- It is likely that the most massive clumps will 
reach M gap , but it will be more difficult to form a clump 
with the mass M[ s . 

From Fig. 4 we see that the most massive clumps can 
reach tens hundreds of solar masses. The expected e-folding 
time for the clump growth is comparable to the orbital pe- 
riod, and thus these high masses will be reached within ~ 20 
orbital periods. The maximum mass would be even larger if 
we included the heating of the disc by external irradiation 
or internal starburst. It is plausible that these very mas- 
sive clumps will form massive stars; the masses of the stars 
may be comparable to the masses of the original clumps; see 
McKee and Tan (2002) and references therein. 

Goodman and Tan (2004) have argued that a proto- 
star accreting from the accretion disc will reach the isolation 
mass; see also Nayakshin and Cuadra (2005) and Nayakshin 
(2005). This conclusion was criticized by Milosavljevic and 
Loeb (ML, 2004), who have shown that the mini-accretion 
disc around the star itself will be unstable to fragmentation. 
Our picture of merging clumps forming a massive protostel- 
lar core is not susceptible to the ML's criticism. 



4.1 Fragmenting disc 

So far we have mostly considered an idealized situation 
where the disc is self-luminous and not irradiated externally. 
However, a gaseous disc in the Galactic Center (when there 
is one) is heated from outside by the light from bright stars 
in the cusp (mostly UV), and re-radiates this in the infrared. 
If the disc is optically thick to its thermal radiation, then 
the temperature T of the disc is uniform in vertical direction 
and is determined by the Stefan-Boltzmann law: 



aT 4 



L-„ 



ttR 2 



(28) 



where R is the disc radius. If the disc is optically thin to its 
thermal radiation, then determining the vertical tempera- 
ture profile of the disc requires more care. The superheated 
dust layer forms on its skin and the interior is heated by the 
infrared light emitted by this layer (Chiang and Goldreich, 
1997). Two cases can be distinguished: first, in which the 
disc interior is thin to its own thermal radiation but thick 
to the one from superheated dust, and second, in which the 
disc is thin to both its own radiation and that of the dust. 
In both cases, however, the mid-plane temperature will be 
slightly, but not much, greater then the estimate in Eq. (28) 
(i.e., by a factor of order r _1//4 in the first case and (es/ti) 1 ^ 4 
in the second case; here r is the disc's optical depth to its 
thermal radiation and e s ,i are the dust emissivities in the 
superheated layer and mid-plane, respectively). We will see 
that the actual optical depth of the disc to its thermal radi- 
ation will be of order 1. 

The current total luminosity of stars is estimated to be 
Lo = 2 x 10 7 Lq (Paumard et al., 2005). Lets assume that a 
fraction (3 of the stellar light is intercepted by the disc (this 
fraction is of order 1 since the disc is optically thick for the 
UV). From Eq. (28), the disc temperature is§ 



(27) T ~ 84/3 1/4 (L/L ) 1/4 (0.5pc/i?) 1/2 K. 



(29) 



and the speed of sound is 

c s ~ 4.5 x 10 4 /3 1/s (L/Lo) 1/8 (0.5pc/i?) 1/4 cm/sec (30) 

From Figure 1 we see that for orbital periods greater 
than 600 years (this corresponds to r > 0.06pc for SgrA*) 
the temperature of self-luminous disc is smaller than that of 
the gaseous disc in the Galactic Center, and thus the disc 
is supported by external irradiation prior to fragmentation. 
At the fragmentation boundary, the Toomre Q parameter 
equals 1: 



Q 



(31) 



a The disc will be also heated by the stars from the central star 
cluster when they collide with the gas (Syer et al., 1981). The rate 
of heating per area can be estimated as ~ Er 2 n<r 3 ~ 5Lq /pc 2 , 
where S is the disc column density, r* is the typical stellar radius, 
and n and a are the density of stars and the velocity dispersion 
of the cluster, respectively. In the numerical estimate we have 
taken the values appropriate for the Galactic Center environment: 
r = r-0, E = 10g/cm 2 , n = 10 6 pc- 3 , and a = lOOOkm/s. Clearly, 
the power input from this process is negligibly small compared to 
that from the external irradiation. 
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and thus from Eq. (30), 

S = 2.7 x (0.5pc/r)- 3/2 (/3L/L ) 1/8 (0.5pc/i?) 1/4 g/cm 2 (32) 

Thus the mass of the fragmenting disc is 

M disc = 4n(R 2 * E(R) - Rf n E(R in )) 

= 2.5 x 10 4 [(i?/0.5pc) 1/2 - (i? in /0.06pc) 1/2 ] x(33) 

(aL/L o ) 1/8 (O.5pc/i?) 1/4 M , (34) 

where Ri n is the inner radius of the disc. I have taken this 
inner radius to be Ri n — 0.06pc, since inside this radius the 
critical column density necessary to achieve fragmentation is 
much larger than outside, see Fig. 2. The value of i?i n agrees 
well with the observed inner edge of the clockwise stellar 
disc, which is 0.03pc in projection on the sky (Paumard et 
al., 2005). The disc mass obtained in the Equation above 
appears to be twice as large as the value estimated from 
observations by Paumard et al. (2005); see also Nayakshin 
et al. (2005). The difference can be explained if either (a) 
the stellar luminosity was smaller by ~ 300 than the current 
one during the disc formation; this is unlikely since even 
though most of the current luminosity is due to the recently- 
formed stars, there is substantial UV flux from post-AGB 
stars in older populations and the center of the galaxy has 
significant soft-x-ray sources which should predate the most 
recent starburst, (b) the starburst occupied the outer half of 
the disc, or (c) only half of the gas was converted into stars, 
and the rest left the GC as wind. 

Equation (32) makes a prediction that £ oc r -1 ' 5 dur- 
ing fragmentation. But what about the number density of 
formed stars? If the typical stellar mass scales with the iso- 
lation mass of the disc (Goodman and Tan 2004, section 3 
of this paper), then 

M isolation oc (£r 2 ) 3/2 oc r ' 75 , (35) 

and the number of stars per unit are scales as 

dN/d(Area) oc r~ 2 ' 25 . (36) 

This is in good agreement with the data in Paumard et 
al. (2005), who find dN/d(Area) oc r~ 2 

4.2 S-stars 

It is not easy to visualize how eccentric orbits orbits of the 
S-stars could be consistent with their birth in a disc, and so 
far the ideas for their origin have invoked stellar dynamics. 
Recently, two interesting suggestions have been put forward. 
Alexander and Livio (2004) have proposed that the S-stars 
got captured by an exchange interaction with stellar-mass 
black holes in the SgrA* cusp; however this scenario seems 
to be disfavored by current observations (Paumard et al, 
2005). Gould and Quillen (2003) have argued that in the 
past the S-stars were members of binaries. When the bi- 
nary passes close to SgrA*, one of its members gets ejected 
and the other remains bound to SgrA* on an eccentric or- 
bit. However, (a) it is unclear how a young binary would 
get onto a plunging orbit during its lifetime, and (b) some 
of the S-stars' eccentricities are not high; for example SI 
and S13 have eccentricities of around 0.4 (Eisenhauer et al., 
2005). It is certainly impossible to capture stars on orbits 
with such low eccentricity by disruption of stellar binaries, 
since the typical eccentricity scales as (1 — e) ~ (m/M) 1 / 3 



r, pc 



Figure 4. The isolation and gap opening masses plotted as a 
function of radius for the critically fragmenting disc without ex- 
ternal sources of heating. The black hole mass is taken to be 
3 x 1O 6 M . 



where m and M is the mass of the binary and SgrA*, re- 
spectively. Therefore, one has to appeal to a fast relaxation 
twice: the first time, to get the binary on a plunging orbit 
when it is far from SgrA*, and the second time, to get the 
star's eccentricity relaxed to a lower value when it is close 
to SgrA*. 

Rauch and Tremaine (1996) have identified a fast relax- 
ation process near supermassive black holes, and called it the 
"resonant relaxation" (RR) . Here we argue that it is poten- 
tially important for the S-stars; this argument has also been 
made independently by Hopman and Alexander (2006)^ 
The characteristic timescale for the eccentricity change due 
to RR is given by 



t RR ~ (0.5/&) 2 



M \ /16M G 
4 x 10 6 Mo ) V m s 



■) (Porb/yr) 



x!0 6 yr,(37) 



where m s is the typical mass of a star in the cusp [which 
is thought to be dominated by stellar-mass black holes in 
the central O.lpc, see Miralda-Escude and Gould (2000)], 
Porb is the orbital period of the star, and f3 s is the nu- 
merical coefficient, found to be approximately 0.5 by Rauch 
and Tremaine. Equation (37) is valid so long as relativis- 
tic precession is longer than the Newtonian precession due 
to the cusp. Otherwise, the expression for tRR has to be 
multiplied by the ratio of the two precession frequencies. 
The longest measured period of the S-star is that of SI, 



^ We acknowledge Scott Tremaine for attracting our attention 
to the RR as a potentially important relaxation mechanism for 
the S-star relaxation, in 2004. 
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which has P or b = 94yr (Eisenhauer et al., 2005, see also 
Ghez et al., 2005). Thus the longest RR timescale for the 
S-star is of order lOOMyr, well within a possible lifetime of 
the main sequence B stars. The RR relaxation timescales are 
shorter for all the other stars, the shortest one being for S2: 
tRR ~ 40Myr. Given that we do not know the typical mass 
of a stellar-mass black hole in the cuspH, we conclude that 
RR is a viable mechanism for relaxation of orbits of 
all of the S-stars within their lifetime. 

However, there is no mechanism known to us which 
would place the young stars on plunging orbits when they 
are further out from the SgrA*, as required in both Gould 
and Quillen (2003) and Alexander and Livio (2004) scenar- 
ios. We therefore suggest that the S-stars were born dur- 
ing the previous starburst or several starbursts, and then 
migrated towards smaller radii through their parent discs. 
This possibility was already pointed out in LB; it is known 
as type-I migration in planetary dynamics. The migration 
velocity through the disc is given by 



2/3 



(38) 



cf. Eqs. (9), (B4) and (B5) of Rafikov (2002), and Ward 
(1986). Here /3 is a numerical factor of order 5 for Q ^> 1, 
and can be significantly larger for Q ~ 1. From Eq. (30) we 
see that the migration timescale for an S-star embedded in 
the disc at O.lpc is of order 10 J yr, i.e. much shorter than the 
star's lifetime. Complication arises if the star opens a gap 
in the disc. However, in this case a fast runaway migration 
is possible due to co-rotation torques acting on the star; the 
speed of the runaway migration is comparable to that of 
type-I migration (Artymowicz 2003, Masset and Papaloizou 
2003). 



If on the other hand, the black hole is not massive enough 
to open the gap, it will migrate inwards by exciting density 
waves in the disc (type-I migration, see previous section). 
The speed of this inward drift is given by Eq. (38). For a 
self-gravitating disc with Q ~ 1, we find 



^inspiral ~ lQ 7 yr 



5 30h 10M, 



o 



P r M b h V O.lpc 



M 7 U 



(40) 



What determines whether the gap is open or not is 
whether a stellar-mass black hole has time to accrete a few 
thousand solar masses of gas, which would put it above 
the gap-opening threshold for the typical disc parameters. 
The Eddington-limited accretion occurs on a timescale of 
~ 10 8 yr, much longer than the characteristic type-I inspiral 
time. Thus, in this case the black holes don't accrete much 
on their way in. On the other hand, the Bondi-Hoyle for- 
mula predicts the mass e-folding timescale of a few hundred 
years. Thus, if the black hole is allowed to accrete at the 
Bondi-Hoyle rate, its mass increases rapidly until it opens 
a gap in the disc. Then, the inspiral proceeds via type-II 
migration. 

From Eqs. (39) and (40) we see that the embedded 
black hole experiencing type-I or type-II migration would 
merge with the central black hole on the timescale 10 6 — 10 7 
years, shorter than the typical timescale of an AGN activity. 
Thus it is plausible that the daughter disc-born black hole is 
brought towards the parent central black hole; the mass of 
the daughter black hole might grow significantly on the way 
in. Gravitational radiation will eventually become the dom- 
inant mechanism driving the inspiral, and the final merger 
will produce copious amount of gravitational waves. In the 
next subsection we show that these waves are detectable by 
LISA for a broad range of the black hole masses and the disc 
accretion rate. 



5 MERGER OF THE CENTRAL BLACK HOLE 
AND THE DISC-BORN BLACK HOLE. 

Stars with masses of a few tens of solar masses will pro- 
duce black holes as the end product of their rapid (~ 10 6 yr) 
evolution; the characteristic mass of these black holes is be- 
lieved to be around WMq. A recent work by Mirabel and 
Rodrigues (2003) shows that the black holes whose progen- 
itors have masses > 40M Q do not receive a velocity kick at 
their birth. Thus, significant fraction of the disc-born black 
holes will remain embedded in the disc. 

It is likely that the newly-born black hole in the disc 
inspirals towards the central black hole. We imagine that the 
stellar-mass black hole is embedded into a massive accretion 
disc which forms due to continuing infall of gas from the 
galactic bulge, after the black hole is born. If the black hole 
opens a gap in the disc, it will move towards the central black 
hole together with the disc (type-II migration; Gould and 
Rix 2000, and Armitage and Natarajan 2001). The timescale 
for such inspiral is the accretion timescale, 

m« M 7~ 1/2 ^0-lpcy 3 / 2 ( r \ 2 
tlMplr(a ~10yr— -(-— ) (—) . (39) 



II Top-heavy IMF may imply top-heavy black holes! 



5.1 Influence of the accretion disc on the inspiral 
signal as seen by LISA 

It is realistic to expect that LISA would follow the last year 
of the inspiral of the disc-born black hole into the central 
black hole. Generally, one must develop a set of templates 
which densely span the parameter space of possible inspiral 
signals. In order for the final inspiral to be detectable, one 
of the templates must follow the signal with the phase-shift 
between the two not exceeding a fraction of a cycle. There- 
fore, if the drag from the disc will alter the waveform by a 
fraction of a cycle over the signal integration time (e.g., 1 
year), detection of the signal with high signal-to- noise ratio 
will become problematic. Below we address the influence of 
the accretion disc on the final inpiral waveform. 

The issue of gas-drag influence on the LISA signal was 
first addressed by Narayan (2000); see also Chakrabarti 
(1996). Narayan's analysis is directly applicable to low- 
luminosity non-radiative quasi-spherical accretion flows, 
which might exist around supermassive black holes when 
the accretion rate is < 0.01 of the Eddington limit. Narayan 
concluded that non-radiative flows will not have any ob- 
servable influence on the gravitational-wave signals seen by 
LISA. Below we extend Narayan's analysis to the case of 
radiative disc-like flows with high accretion rate, which are 
likely to be present in high-luminosity AGNs. 
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Consider a non-rotating central black hole of mass 
Ma x 1O 6 M0, accreting at a significant fraction ra of the Ed- 
dington rate, M — mM cdd . The accretion disc in the region 
of interest (< WRs, where R s is the Schwartzschild radius of 
the central black hole) is radiation-pressure dominated, and 
the opacity k = 0.4cm 2 /g is due to the Thompson scatter- 



ing. By following the standard thin-disc theory 
and Sunyaev, 1973), we get 



(Shakura 



64ttc 2 



3/2 



S(r) = """" . ~ 4g/cm 2 -V ( — ) ' , (41) 

where e is the efficiency with which the accreted mass con- 
verts into radiation, and r s is the Schwartzschild radius of 
the central black hole; and 



3Mk 
8ttc 



(42) 



The disc black hole in orbit around the central black 
hole excites density waves in the disc (Goldrcich and 
Tremaine, 1980); these waves carry angular momentum flux 



F ~ (GM bh ) 2 



(43) 



here, as before, M b h is the mass of the orbiting disc black 
hole. Ward (1987) has argued that the torque acting on the 
orbiting body is dL dw /dt ~ (h/r)Fo. We can compute the 
characteristic timescale for the orbit evolution due to the 
density-waves torque: 



1 M M 



dL d „/dt fi M bh Sr 2 



(44) 



One must compare idw to the timescale t gw of orbital evo- 
lution due to gravitational-radiation-reaction torque: 



8i„ 



5 or. 



8 GM b 



(0 



(45) 



Here, t m is the time left before the disc and central holes 
merge, i.e. the integration time for the LISA signal. Opti- 
mistically we could expect to follow the LISA signal to 0.1 
of a cycle. Therefore, if q = 10n m i gw /idw is less than unity, 
the disc drag does not impact detection of the final inspiral; 
cf. Eq. (16) of Narayan (2000). Here n m = f2i gw /(57r) is the 
number of cycles the disc black hole will make before merg- 
ing with the central black hole. Using Eqs. (41), (42), (44), 
and (45), we get 



q ~ 2 x 10 



4, 



(Mbh/IOM©) 13 / 8 



.21/8 



M, 



13/4 



(46) 



where e = O.leo.i, a = O.Iqo.i, and t m — lyr*i mi!/r . We see 
that for a large range of parameters q < 1, and the disc drag 
does not influence the inpiral signal. However, note that q 
is a very sensitive function of m and Mbh. For instance, the 
disc will influence significantly an inspiraling 100Mq black 
hole when the accretion rate is down to a few percent of 
the Eddington limit. One then needs to reduce the influence 
of the disc on the LISA signal by choosing to observe the 
smaller portion of the final inspiral, i.e. by choosing a smaller 
integration time t m . 

** The disc is no longer thin close to the central black hole, 
however our estimates of the disc structure should be correct to 
an order of magnitude. 



There is another important source of drag experienced 
by the inspiraling black hole; it was first analyzed by 
Chakrabarti (1993, 1996). Generally, there is a radial pres- 
sure gradient in an accretion disc; this pressure gradient 
makes the azimuthal velocity of the disc gas slightly dif- 
ferent from a velocity of the test particle on a circular orbit 
at the same radius. Therefore the inspiraling black hole will 
experience a head wind from a gas in the accretion disc; by 
accreting gas from the disc the black hole will experience the 
braking torque which will make it lose it's specific angular 
momentum. This torque is given by 



Twind = M bh Avr, 



(47) 



where Mbh is rate of accretion onto the inspiraling black hole 
from the disc, and Av ~ c 2 s /v is the speed of the headwind 
experienced by the black hole moving with the orbital speed 
v . We assume that the inspiraling hole accretes with the 
Bondi-Hoyle rate, 



np(GM bh f/c 3 s , 



(48) 



where p = E//i = £f2/c s is the density of the ambient disc 
gas. By using the last equation in Eq. (47), we obtain the 
expression for the torque r w i n d, which turns out to be the 
same as the torque from the density waves, up to the nu- 
merical factor between 1 and 10. We see therefore that in- 
clusion of Chakrabarti's "accretion" torque is important for 
the detailed analysis, but does not qualitatively change our 
conclusions. 

Can the orbiting black hole open a gap in the disc during 
it's final inspiral? In order to overcome the viscous stresses 
which oppose opening the gap, the mass of the orbiting black 
hole should be greater than the threshold value given by 

Mbhgap =± ^/4M>/0 2 ' 5 M, (49) 

sec Eqs. (4) of Lin and Papaloizou (1986). The radius n n 
from which the inspiral begins is given by 



r in = 4(M b h/lOM ) 1/4 M 6 r 1/2 4 / ,t /r .r s 



(50) 



and the scaleheight of the disc is radius-independent in the 
radiation-dominated inner region (this is true only if you 
treat the central black hole as a Newtonian object): 

h=^r s . (51) 

The gap will be open in the disc, therefore, if the mass of 
the inspiraling black hole exceeds 



^20/13 M -2/13 lo 4 M0 



(52) 



During the final year of merger of the IOMq disc-born black 
hole and the 10 6 M o( jot central black hole, the gap will be 
open if the accretion rate onto the central black hole is 
about five percent of the Eddington limit. Such gap-opening 
merger can result in sudden changes in the AGN luminos- 
ity and produce an optical counterpart to the gravitational- 
wave signal. 

5.2 The merger event rate as seen by LISA 

The details of black-hole formation and evolution in the 
disc are uncertain, and it seems impossible to make a re- 
liable estimate of an event rate for LISA from this chan- 
nel of black-hole mergers. Nonetheless, the expected top- 
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heavy IMF gives us reason for optimism. It is worth work- 
ing through a simple example to show that the merger of 
disc-born holes with the central hole may be an important 
source for LISA. 

Studies of integrated light coming from Galactic Nuclei 
show that the supermassive black holes acquire significant 
part and perhaps almost all of their mass via accretion of 
gas; see, e.g., Yu and Tremaine (2002) and references therein. 
Assume, for the sake of our example, that a mass fraction 
77 of this gas is converted into lOOM© black holes on the 
way in, and that all of these black holes eventually merge 
with the central black hole. LISA can detect such mergers 
to z = 1 if the central black hole is between 10 5 and 10 7 
solar masses (one needs to assume that the central black 
hole is rapidly rotating at the high-mass end of this range). 
The mass density of such black holes in the local universe 
is ~ 10 5 M©/Mpc 3 (Salucci et. al. 1999). We can estimate, 
using Figure 3 of Pei et. al. (1995), that about 10 percent of 
integrated radiation from Galactic Nuclei comes from red- 
shifts accessible to LISA, z < 1. This implies that supermas- 
sive black holes acquired 10 percent of their mass at z < 1. 
We shall therefore assume^ that 10 percent of the mass of 
the black holes in LISA mass range was accreted at z < 1. 
This implies that there were ~ IOO77 = 770.01 mergers per 
mpc 3 which are potentially detectable by LISA. When we 
multiply this by volume out to z = 1 and divide it by the 
Hubble time, we get an estimate of the LISA event rate from 
such mergers, 

dN/dt toyulodcl ~ 10770.01/yr. (53) 

The estimate above should be treated as an illustration of 
importance for our channel of the mergers, rather than as a 
concrete prediction for LISA. 

Currently, it is not known how to compute well the grav- 
itational waveform produced by an inspiral with an arbitrary 
eccentricity and inclination relative to the central black hole. 
This might pose a great challenge to the LISA data analysis. 
Indeed, in the currently popular astrophysical scenario, the 
stellar-mass black hole gets captured on a highly eccentric 
orbit with arbitrary inclination relative to the central black 
hole (Sigurdsson and Rees, 1997). By contrast, in our sce- 
nario the inspiral occurs in the equatorial plane, and the sig- 
nal is well understood (Hughes 2001 and references therein). 
Circular inspiral is also expected in the binary-disruption 
scenario of Miller et al. 2005; however there the inspiral is 
not expected to be preferentially aligned with the equato- 
rial plane. The template for detection of a circular inspiral 
is readily available, and the merger channel we consider has 
a clear observational signature which distinguishes it from 
other channels. It is an open question whether the inspiral- 
ing hole can acquire high eccentricity by interacting with 
the disc or with other orbiting masses (Goldreich and Sari 
2003, Chiang et. al. 2002); it seems likely that at least in 
some cases interaction with the disc will act to circularize 



the orbit. Finn and Thorne (2000) have performed a detailed 
census of parameter space for circular-orbit equatorial inspi- 
rals as seen by LISA. 



6 CONCLUSIONS 

In this paper, we rely on numerical simulations by Gammie 
(2001) to develop a formalism for self-gravitating thin discs 
which are gaining mass by continuous infall. We present a 
way to calculate the critical temperature, surface density, 
and scaleheigh of the disc just prior to fragmentation. Our 
formalism naturally includes both optically thin and opti- 
cally thick discs. 

We then speculate on the outcome of the nonlinear 
physical processes which follow fragmentation: accretion and 
merger of smaller fragments into bigger ones. We find an up- 
per bound on the mass of final, merged fragments, and we 
give a plausibility argument that some fragments will indeed 
reach this upper bound. We thus predict that very massive 
stars of tens or even hundreds of solar masses will be pro- 
duced in self-gravitating discs around supermassive black 
holes. This picture explains well recent observations of the 
stellar discs in the Galactic Center, and we clarify how the 
peculiar orbits of the S-stars may be generated by a com- 
bination of inward migration through the gaseous disc and 
Resonant Relaxation in the central arc-second. 

The end product of fast evolution of the massive stars 
will be stellar-mass black holes. We make an argument that 
in the case when accretion is prolonged, the disc-born black 
holes in AGNs find a way to merge with central black holes. 
We consider a purely toy-model example of what the rate 
of such mergers might be, as seen by LISA; we illustrate 
this rate might be high enough to be interesting for future 
gravitational-wave (GW) observations. The GW signal from 
this merger channel is distinct from that of other channels, 
and can be readily modeled using our current theoretical 
understanding of the final stages of the inspiral driven by 
gravitational-radiation reaction. We show that for a broad 
range of accretion rates and black-hole masses the drag from 
accretion disc will not be large enough to pollute the signal 
and make inspiral template invalid. In some cases, the inspi- 
raling hole will open a gap in the accretion disc close to the 
central black hole, thus producing a possible optical coun- 
terpart for the gravitational-wave burst generated by the 
merger. 
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tt This assumption has an extra caveat for the low-luminosity 
AGNs, significant fraction of which may be powered by tidal dis- 
ruption of stars (Eracleous et al. 1995, Milosavljevic et al. 2006). 
We note that in one low-luminosity AGN, in the nucleus of Circi- 
nus galaxy, there is a firm evidence of an extended accretion disc 
(Greenhill et al. 2003). 
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